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ABSTRACT 

In chameleon field theories a scalar field can couple to matter with gravitational strength and still 
evade local gravity constraints due to a combination of self-interactions and the couplings to matter. 
Originally, these theories were proposed with a constant coupling to matter, however, the chameleon 
mechanism also extends to the case where the coupling becomes field-dependent. We study the 
cosmology of chameleon models with power-law couplings and power-law potentials. It is found that 
these generalized chameleons, when viable, have a background expansion very close to ACDM, but 
can in some special cases enhance the growth of the linear perturbations at low redshifts. For the 
models we consider it is found that this region of the parameter space is ruled out by local gravity 
constraints. Imposing a coupling to dark matter only, the local constraints are avoided, and it is 
possible to have observable signatures on the linear matter perturbations. 

Subject headings: (cosmology:) cosmic background radiation; cosmology: miscellaneous; cosmology: 
observations; cosmology: theory; cosmology: large-scale structure of universe 



1. INTRODUCTION 

The origin of dark energy (DE) responsible for the cos- 
mic acceleration remains a mystery. A host of indepen- 
dent observations have supported the existence of dark 
energy over the past decade, however, no strong evidence 
was found yet implying that dynamical DE models are 
better than a cosmological constant (see e.g. IDavis et alj 
20071) for results on observation of dark ene rgy and 



Copeland et~alll2006l : iDurrer fc Maartensll2008D for the- 



oretical overviews) . The first step towards understanding 
the origin of DE would be to detect some clear deviation 
from the ACDM model observationally and experimen- 
tally. 

Scalar fields (quintessence) are natural DE candi- 
dates, with an equatio n of state different from — 1. 
In qu intessence models (jWang et al.l 120001 : iZlatev et al.l 
Il999f) the scalar field is slowly rolling down its potential, 
its energy density is dominated by the potential energy 
and almost remaining constant provided that the poten- 
tial is flat enough. However, this means that the mass 
of the scalar field is in general very light and as a result 
the scalar field almost does not cluster so that its effects 
in cosmology are mainly on the (modified) background 
expansion rate. 

However, a certain class of theories have been pro- 
posed, in which the scalar field(s) properties depend 
on the environment: these are the class of chameleo n 
field theories, proposed in iKhourv fc Weltmanl (|2004l ). 
that employed a combination of self-interactions of the 
scalar-field and couplings to matter to avoid the most 
restrictive of the current bounds. In the models that 
they proposed, a scalar field couples to matter with 
gravitational strength, in harmony with general expec- 
tations from string theory, whilst, at the same time, re- 
maining relatively light on cosmological scales. Origi- 
nally, the chameleon was proposed with a co nstant cou- 
pling t o matter, however, it was shown in IBrax et al.l 
(|2010bD that the chameleon mechanism extends to the 
case where the coupling becomes field-dependent. The 



chameleon mechanism ca n also be present i n other types 
of models like shown in IWei fe Call (|2005f ) for the case 
of k-essence. Because chameleons are allowed to cou- 
ple strongly to matter and at the same time exhibit a 
long ranged fifth-force in space, they could leave a strong 
signature on the growth of the matter perturbations. 
Chameleon models also have signatures in local gravity 
experiment s which can be within reach of near-future ex- 
periments (IKhourv fc Weltmanl 120041: IBrax et al.ll2010al : 
ISteffen fc Gammev Collaborationll200cfh 

The modified evolution of the matter density pertur- 
bations can provide an important tool to distinguish 
generally modified gravity DE mode ls, from DE models 
inside GR like the A C DM model dCarroll et al.1 |200F 
Faulkner etHI [20071: ISong et all l2007at 



20071: ISong et all I2007bl: IPogpsian fc Silvestrl pi)8 
Tatekawa fc Tsuiikawal l2008t iMota fc Barrowf 12004 




Ovaizu et all l2008t iKovama et al.l 120091: IMota et al 



2007 : iKoivisto fc Motall2008t iBourliot et al 1 120071: IMota! 



2008D . In fact the effective gravitational constant 
G c ff which appears in the source term driving the 
evolution of matter perturbations can change signif- 
icantly relative to the gravitational constant G. A 
useful way to describe the perturbations is to write 



the growth function / 



d log <5„ 
log a 



as / = f2 m (z) 7 where 



Q m is the density paramete r of non-relativistic matter 
(bary onic and dark matter) (jPeebles! 119841 : ILahav et al.l 
I1991D . One has 7 ~ 0.55 at low redshifts in the 
ACDM-model dWang fc Steinhardtl [19981 iLinderl [20051 : 
iHuterer fc Linderl I2007D . While 7 is quasi-constant in 
standard (non-interacting) DE models inside GR, this 
needs not be the case in modified gravity models. An 
additional important point is whether 7 can exhibit 
scale dependence (dispersion). When this happens the 
resulting matter power spectrum acquires an additional 
scale dependence wh i ch is n ot found in ACDM. 

In iGannouii et al.l (|2010f ) a chameleon model, with a 
constant coupling, was constructed where the present 
value of the growth index 7 can be as small as 7 = 0.2 
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together with a significant redshift dependence. This al- 
lows to clearly discriminate this model from ACDM. 

In this paper, we investigate the cosmological proper- 
ties of chameleon models where the matter-coupling is 
a power-law function of the scalar-field. This article is 
divided in 5 main sections. In section 2 we give a short 
review of chameleon models and define the models we will 
be looking more closely at. In section 3 we discuss the 
background evolution. In section 4 we discuss local and 
cosmological constraints on our models, and in section 5 
we consider the evolution of the matter perturbations in 
our models. 

2. CHAMELEON THEORIES 

In this section we review the basic properties of 
chameleon models. We start by considering the scalar- 
tensor theory described by the action 



S 



where M p \ 



RM, 



1 



A) 
(1) 



is the reduced Planck scale, 
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the determinant of the metric g^, R is the Ricci-scalar 
with respect to g^ and ipi are the different matter-fields 
(including radiation). The matter field i couples to the 

~(i) 

Jordan-frame metric g^l which is related to the Einstein- 
frame metric g^ via a conformal rescaling on the form 



g(i) 
.y [iv 



C2) 



In the following we focus, for simplicity, on the case where 
all the matter-fields couple to <f> with the same (3i(<f>) = 
f3((f>). Thus, a general model of this kind is then uniquely 
specified by stating the two functions V(<j>) and (3(4>). 

A variation of the action Eq. ((T|) with respect to (j> 
yields the field-equation 



n<j) 



where T^ v 



; ac, 



(3) 



- is the energy-momentum ten- 
v -g a 9v b 
sor of the matter fields. The energy-momentum ten- 
sor for radiation vanishes meaning that chameleons do 
not couple directly to photons. However, a coupling to 
photo ns, which have some interesting observable signa- 
tures (IDavis et al.ll2009l : ISchelpdl20la IDavis et al.ll2010l : 
iBrax fe Zioutasl 120101) . can be introduced by modifying 
the electromagnetic field strength F% v — > e^ 6 'F^ v . 

In the perfect fluid approximation, we have g^JT^ = 
—p m where p m is the Jordan-frame energy density of 
non-relativistic matter. The density p m is conserved 
with respect to the Jordan- frame metric g^ v . In the 
Einstein- frame the density p rn = /5 m e 3/3 ^- ) is both con- 
served and ^-independent. The field equation, in the 
Einstein-frame, can be written 



□0 = V c s,6 



(4) 



where V e g is the effective potential. In the following, the 
quantity will be referred to as the coupling and P((f>) 
as the coupling-function. We will also use the notation 

44 and A, 6 =i A 



2.1. Chameleon thin-shell mechanism 

We will in this section give a short review of the 
chameleon thin-shell mechanism in which these theories 
can have a strong matter coupling and still be viable. 

The coupling Eq. © of <fi to matter leads to the in- 
troduction of a fifth-force in nature. In linear theo- 
ries of massive scalar fields the superposition principle 
holds, meaning that the larger a massive body is the 
stronger the fifth-force becomes. Chameleon theories be- 
have quite opposite: In situations where massive bodies 
are involved, the chameleon field is trapped inside the 
bodies and its influence on other bodies is only due to a 
thin shell close to the surface of the bodies. This leads 
to a shielded fifth- force which becomes hard to detect. 

In order to have a chameleon mechanism in a theory 
the following conditions must be satisfied 

• For a given density po the effective potential V e s = 
V{4>) + p e PW must exhibit a minimum <p m \„. 



• The mass, m m i n = \/V e g , c fxt>(4 > min) , of small oscilla- 
tions about this minimum must be a real increasing 
function of pq. 

and puts the following constraints on the potential V((j)) 
and coupling-function /?(</>): 



i) V„j> (3,4, < 0, ii) V,<j, > 0, 



iii) /3,00<* f3,<t, > 
(5) 

Note that the second condition above is required if we 
want the chameleon to exhibit a long-ranged force in the 
vacuum of space. The final requirement, which is the 
most important, is that there exists a chameleon thin- 
shell mechanism in which the fifth-force, for sufficiently 
large bodies, becomes suppressed relative to the gravita- 
tional force. 

To see how this mechanism works in practice, we con- 
sider a spherical object with constant density p Cl radii R c 
and mass M c embedded in a background of homogeneous 
density p^ -C p c - In such a static spherical symmetric 
space time, the field equation reads 



dr 2 



air 



V,6 +p c (3,6 e 



0M 



r<R c 
r > R r 



(6) 



The effective potential have two minimums at <f> = <f> c 
and 4> = 4>b satisfying 

V l6 (4> c ) + fa W> c )p c e«« = (7) 
V, 6 (&) + p, 6 (<h)Pbe^ = (8) 
The (^-mediated fifth-force per unit mass is given by 

F 6 = -p, 6 V<j> for r>R c (9) 

In the following we will consider the simplest case /3(</>) = 
■g^-, in which the coupling P,6= Q/M p \ is constant, 
together with an arbitrary runaway potential like e.g. 

V(*) = 

For this case, the field profile outside t he body have 
been derived in iKhourv fc Weltmanl (|2004f ) and reads 



<2cff M c 



pi 



(10) 
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where = ^JVoS^J^M IS the mass of the field in the 
background and Q c ff is the effective coupling. The value 
for Q ff depends on the value of the thin-shell factor 



eth 



4>b 



6QA/ P i$ c 



(11) 



where $ c is the Newtonian gravitational potential. When 
the body is small in the sense that eth > 1 the field acts 
like a linear scalar field and Q c r = Q. 

For large bodies, e t h <C 1, the field will be stuck at the 
minimum <j) c inside the body, and the only changes in <j> 
takes place in a thin-shell close to the surface. We say 
that the body has developed a thin-shell. 

In this case we find Q c r = 3Qe t h- Thus, the cou- 
pling strength is suppressed relative to Q. More accu- 
rate formulas for Q f ff whe n e t h < 1 can be found in 
iTamaki fe Tsuiikawal (|2008Q , 

The amplitude of the fifth-force on a test-particle of 
unit mass outside the body is 



F,/, — 2QQ c tf 



GAL 



for r < m. 



(12) 



For a sufficiently large body we have 2QQ c r <C 1, and 
the resulting fifth-force is suppressed relative to the grav- 
itational force. When test-masses used in local gravity 
experiments have thin-shells, the experimental bounds 
are easily evaded. 

For more general coupling-functions /?(</>) the analysis 
above can be much more complicated due to the cou- 
pling ,6 in Eq. ((9J being field-dependent, see lBrax et al.1 
(|2010bf )). We present a derivation of the thin-shell solu- 
tion for the power-law coupling /3((f>) z 



m > 1 in section 4. 



with 



2.2. Our models 

We will look more closely on the following models, in 
which the cosmological constant is part of the potential: 



Model A 



V((/>) = M 4 

m = 



i 



- I M. 



f) 



(13) 



with n > and m > 1. This model is a generalization 
of the original chameleon model (m = 1) introduced in 
IKhourv fc Weltmanl (12001 ). 



Model B 



V(<(>) = M 4 



i 



(*)" 



(14) 



with n > 2, m > and Hq being the current Hubble 
parameter. This model was introduced, and a wide range 
of local gravity constraints was calculated, in lBrax et alJ 
(|20T0bTl . 

In order for the chameleon to act as dark energy we 
need to choose M 4 ~ A (M ~ 10 _3 eV). Note that this 
value is also required in order to satisfy local gravity 
constraints. Thus, the fine-tuning problem of the cosmo- 
logical constant is also present in these models. 

3. BACKGROUND EVOLUTION 

In this section we discuss the background cosmological 
evolution of chameleon models. 



We consider a flat Friedmann-Lemaitare-Robertson- 
Walker (FLRW) background metric 

ds 2 = -dt 2 + a(tf{dx 2 + dy 2 + dz 2 ) (15) 
The corresponding background equations are given by 

3tf 2 A/ p 2 ! = p m e^ + Pr + X -4> 2 + V{4>) (16) 
p m + 3H Pm = (17) 

pr + AH Pr = (18) 

The field equation for <j> becomes 



+ 3H(f> + V cS ^ = 
V eS = V{4>) + Pm e^ 

We also introduce the density parameters 



(19) 



- - 1 , 



Pr 

beg 
- K±j£ 



(20) 



The background ev olution of the original ch ameleon 
model introduced in IKhourv fc Weltmanl (|2004f) , which 
corresponds to Mod el A Eq. (|T3")) with m = 1, was thor- 
oughly discussed in lBrax et al.l (|2004bf ). We will in the 
next section show that the chameleon behaves very sim- 
ilar in the general setting. There are however some im- 
portant differences. 

3.1. Attractor solution 

We show the existence of an attractor solution where 
the chameleon follow the minimum of its effective poten- 
tial </> = </> m i n (i) as long as the condition ^> H is 
satisfied. 

Suppose the field is at the minimum at some time ti. 
Then a time later due to the red shifting of the mat- 
ter density the minimum 4> m i n has moved to a slightly 
larger value (or smaller, depending on the form of the 
coupling). The characteristic timescale for this evolution 
is the Hubble time 1/H. Meanwhile the characteristic 
timescale of the evolution of <j> is given by When 
<C H the response-time of the chameleon is much 
larger than 1/H, the chameleon cannot follow the min- 
imum and starts to lag behind. But if 3> H then 
the response-time of the chameleon is much smaller than 
1/H, the chameleon will adjust itself and adiabatically 
start to oscillate about the minimum. This can also be 
seen from the analogy of Eq. (fl9|) with a driven harmonic 
oscillator 



2C,ujx + ui"x = 



(21) 



This equation will have a solution which oscillates with 
a decreasing amplitude as long as Q < 1. When the field 
is close to the minimum we can approximate 



3H(f> + mi 



= 



(22) 



The condition £ < 1 reduces to -^j- > 1, and we will 
require that m 2 3> H 2 is satisfied from the early universe 
and until the present era. 
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3.2. Reaching the attractor 

We discuss how the chameleon acts in the early uni- 
verse, and how the field converges to the attractor. For 
simplicity we will only focus on the case where (/> m ; n < 
(Model B Eq. (HU)), the case fa nin > (Model A Eq.[T3J|) 
is analogous. 

If we release the field at some time ti, at some initial 
value fa, we can have two cases: 

3.2.1. Undershooting: (pi S> <^ m in(ti) 
In this case the field equation can be approximated by 



m<p « -v, 4 



(23) 



which is the same equation as in quintessence. The 
driving term will dominate over the friction term when 
3> H 2 and will start driving the field down towards 
the minimum. If this condition is not satisfied, the field 
will be fixed at <pi until the Hubble factor has had time to 
be sufficiently redshifted allowing the field to start rolling 
down the potential. 

When the field starts to roll it drops to, and go past, 
<P — fanin- Here the approximation Eq. (|23|) cannot be 
used anymore, but the field will usually have to much 
kinetic energy to settle at the minimum and will be 
driven past it. The further below the minimum the field 
is driven, the larger the factor (3,^ p m becomes. Even 
though p m is very small in the radiation era it will even- 
tually kick in and drive the field up again. More impor- 
tantly, we will also have a contribution from the decou- 
pling of relativistic matter which will be discussed in the 
next section. This will make the field oscillate around 
the minimum, and as long as m^/H 2 ^> 1 the ampli- 
tude of the oscillations will be damped, making sure that 
the field quickly converges to the attractor. Thi s can 
be showed explicitly as done in iBrax et al.l (|2004bf ). the 
derivation found there is general and applies to our case 
as well. 

3.2.2. Overshooting: pi <C p m in(ti) 

In this case the potential term V,<f, can be ignored and 
the (^-equation becomes 



+ 3i/0«/3, T; 



(24) 



where we have restored the trace of the energy- 
momentum (EM) tensor. In the radiation-dominated 
era this trace is very small since radiation does not con- 
tribute to the trace, and the field will be frozen at its 
initial value. As the universe expands and cools the dif- 
ferent matter-species decouple from the radiation heat 
bath when the temperature is of the same order as the 
mass of the matter-particles. This gives rise to a trace- 
anomaly where the trace of the EM-tensor gets non-zero 
for about one e-fold of expansion leading to a 'kick' in the 
chameleon pushing it to larger field- values. This trace 
can be written for a single matter-species, see IBrax et al.l 
(I2004bn . as 

Tf* - -^J^-^r^/T) (25) 

where is the effective number of relativistic degrees 
of freedom, gi, Ti and rrii is the degrees of freedom the 



temperature and the mass of species i respectively. The 
r-function is given by 



t(x) 



\Jy 2 - x 2 dy 
e» ±1 



(26) 



and ± refers to bosons and fermions respectively. 

See Fig. Q for a plot of T^/{3H 2 M^) in the radiation 
era. The plot shows that each kick contributes to the field 




Log 10 z 



Fig. 1. 



The trace of the EM-tensor, f2„ 



= -Tjt/^M^, 
in the radiation dominated era for all the different matter species 
decoupling from the radiation heat bath. 



equation approximately as an effective matter-density 
f^m off ~ O(0. 01). By using a delta-function source as 
the kick we can show that the result is to push the field 
up a distance 



\A(P\ = O |/3,^ |M. 



9 * (mi] 



(27) 



where fa is the field- value before the kick sets in. When 
= constant, the field will be kicked almost the same 
amount each time a new species freezes out. But when 
7^ constant the smaller the initial value fa the more 
effective the kick is in bringing it back up. When above 
the minimum the resulting kicks will be balanced by the 
term V,$ which drives the field down again making the 
field oscillate above the minimum before eventually set- 
tling down. 

Due to this mechanism, the chameleon can have initial 
values far below the minimum and still be able to get 
close to the attractor relative quickly. The initial value 
will of course depend on the how the chameleon behaves 
under inflation. If the chameleon couples to the infla- 
ton and sits at the minimum at the onset of inflation, 
then after the inflaton decays to reheat the universe the 
density of matter-species coupled to the chameleon will 
decrease rapidly since most of the energy will go to ra- 
diation. This will lead to a release of the field at a value 
well above the minimum, where the undershoot solution 
applies. As long as H , the field will typically settle 

at the minimum before the onset of the Big Bang Nucle- 
osynthesis (BBN). See Fig. ([2]) for a typical evolution of 
p in the early universe with and without the inclusion of 
the kicks. 

3.3. Dynamics of <p> along the attractor 



■5 




Log 10 z 



Fig. 2. — The evolution of <j> in the early universe as a function of 
the redshift z for two types of initial conditions: <f>i 3> <p m i a (above) 
and <f>i -C m in (below). In this example V(4>) = M 4 + </> 4 , fi{<j>) = 
The dashed (full) line corresponds to the solution where we 
neglect (include) the kicks. The kicks-solution does not reach the 
minimum until z rs 10 s , but due to the large mass of the field it 
starts to follow it right away. 



When the field follows the attractor </> as <j> r 
VeS,<f> ~ 0. Taking the time-derivative yields 



-3H 



V, 4 



m 4, 

for f3(4>) < 1. We further find 

2V ~ 2ml r 



15 we have 



(28) 



where 



Vrn% 
VI 



(29) 



(30) 



As long as r > 1, the field will be slow-rolling whenever 
the condition rr& 3> H 2 is satisfied. 

The equation of state for a minimal coupled scalar 
field (quintessence) is given by cj^ = f„T^ ■ Since 



2 +2V 



the chameleon is not minimal coupled, the time evo- 
lution of p<j, must be computed directly from P4>/p<t> = 
-3H(1 + w cff ). Using Eq. (28]), we find 



W e ff 



1 



1 

r 



(31) 



From this equation we see that the chameleon acts as a 



dark energy fluid as long as | <T, and its only for models 
where T ~ O(l) where we can have a significant deviation 
from lu = —1. See Fig. ([3]) for a typical evolution of uj e g. 



CO 



eff 



Log 



10 



Fig. 3. — The effective equation of state for the chameleon when 
V{4>) = M 4 + <j> 4 , /3(<f>) = Za.. We have u = = -0.2 

during the period before BBN. After the transition to (f> < M, 
the chameleon behaves like a cosmological constant with ui = — 1. 
The oscillations in ui a ff comes from the field oscillating around the 
minimum before eventually settling down. 



(32) 



Model A Eq. (JT3J) yields 

to (n + m) ( c/) min \ " 



r = 1 + — 

n 



n 



M 



The transition from 
for a redshift 



(1 + *) 



(MA' 
mflmo \M\J 



) m ; n < M to </> m in > M takes place 

1/3 



10 1 



'A"" 1 / 3 (33) 



which for typical values of (to, A) is before the time of 
BBN (z as 10 9 ). This means that in the background 
today we have <j> 3> M and therefore V(<fi) as M 4 and 

r> i. 

For Model B Eq. (JUJ we find 



„ to in + to) 

r = i + — + ^ i 

n n 



M 



(34) 



The transition from 
for a redshift 



(l + z) 



4>Q 



> M to m i n < M takes place 
M 



10 



10m^-f 



(35) 



where we have used 



%i - 10 30 . For typical values 

of (to, A) we find a redshift which also is before the time 
of BBN. Today this translates into <fi <Si M giving V((j>) as 
M A and r > 1. 

Thus the models we consider here, will have a back- 
ground evolution very close to that of ACDM. 

3.4. Statefinder parameters 

The statefinder diagnostics, introduced in lSahni et al.l 
(|200l . can be a useful tool for distinguishing different 
DE models. The statefinder parameters are defined by 



r - 1 



3(9-1/2) 



(36) 
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where q = — -^i is the deceleration parameter. Upon 
defining h = H 2 . it follows 



h' 



h" 3h' , . 

r = 1+ 2jT2* < 37 > 



where a prime denoted a derivative relative to x = 
— log(a). ACDM, neglecting the contribution from ra- 
diation, corresponds to the fixed point (r, s) = (1,0). 
For a general chameleon (in units of M p i = 1) we find 



+ 2 (< 
3 V 



' \2 



-il m f)^<j) + — 



1 

3 V 



') 2 h 



V 



1 



2V 



(38) 
(39) 



When the chameleon is slow rolling along the minimum 

and ml > H 2 , 

(40) 
(41) 



of the effective potential, that is <\> m 
we can simplify the above equations to 

r - 1 



2 r * 



3 

T 



When r > 1 we find r ~ 1 and s « 0. It is only 
for r ~ 0(1) that we can have an (observable) devi- 
ation from ACDM. For the models considered here, 
T ^ > 1 and the statefinders will be the same as in ACDM. 
See IGannouii et al.l (|2010D for a viable chameleon model 
where the statefinder parameters deviates significantly 
from ACDM. 

4. CONSTRAINTS 
4.1. Local gravity constraints 

Experimental t ests o f gen eral relativity in the solar 
system, see e.g. IWil (fl99l and searches for a fifth- 
force in nature (jKapner et al l 120071 : iHovle et afl l2001h 
gives strong constraints on any new interactions. For 
chameleon models, these constraints are usually avoided 
due to the chameleon mechanism as long as typical test- 
masses used in the experiments have thin-shells. 

Local gravity cons t raints for Model B Eq. (fT4|) was 
found in IBrax et al.l (|2010b[ ). see Fig. (|4|) for the con- 
straints when ttl = 1 and n = 6, 10. 

In order for the chameleon to affect large scale struc- 
ture formation we need the coupling, evaluated in the 
cosmological background Q^q = |/3,0 today M p \\, to be of 
gravitational strength: 0(1) < Q$q. 

For Model B, this corresponds to the r.h.s of the dashed 
line in Fig. (0)). From this figure we see that this requires 
M < 10" 3 eV as claimed below Eq. (fT4| . 

Local gravity constraints for the original chameleon 
model (Model A with m = 1) have been calcu- 
lated in several papers , see e. g. IKhourv fe Weltmanl 
(l200l: IMotafe Shawl (l200l: 
Ste ffen fc Gammev Collaboratior] 
(I2007T): IBrax et all (|2007f ): iCies et all 



H200lT 



Brax et all 
( 200al h 



IMota fc Shawl (j200fit ). 

Below, we derive thin-shell solutions for Model A 
Eq. ([in]), and use these solutions to find the the local 
constraints coming from tests of the equivalence princi- 
ple in the solar-system. 



4.1.1. Thin-shell solutions for the power-law coupling 
We show the existence of thin-shell solutions for the 
power-law coupling-function /?(</>) = (xty) together 
with an arbitrary run-away potential like e.g. V((f>) = 
M ^ (Mod el A). This model was given a treatment in 

IBrax et all (|2004al ) and it was found that there do not 
exist thin-shell solutions for m > 1 when A = 0(1). For 
values A > 1, which is equivalent to a coupling scale 
M* = <C Mpi, there can indeed exist thin-shells. 
Another conclusion the authors of IBrax et all (|2004a| ) 
reached was the existence of a possible singularity in the 
field-profile. We have calculated the field-profile numer- 
ically and no such singularity was found. The solution 
corresponding to the claimed singularity seems to come 
from a mathematical correct, but non-physical, solution 
to the field equation. 

The field equation in a static spherical symmetric met- 
ric with weak gravity reads 



air 2 



2d(j) 
r dr 



V,<f, +p0,<{, 



(42) 



where we have assumed f3(<f>) <C 1. We consider a body 
with constant density p Cl radii R c and mass M c embed- 
ded in a background of homogeneous density pb and im- 
pose the boundary conditions 



dr 



r=0 



dr 

r—oo 

4>(r — > 00) 



(43) 



Lets start by considering a test body. The field inside 
and outside the body is then just a small perturbation in 
the background (fy^. Solving the linearized field equations 
we find 



4>b 



Qb M, 



4ttM p i 

Qb = P,0 b Mpi 



c e —mtr 



r>R r . 



(44) 
(45) 



The amplitude of the fifth-force between two test bodies, 
located in the vacuum of space, is 

F(r) = 2Q^F gravity (r) for m b r < 1 (46) 

Now lets see what happens for mac roscopic bodies. The 
thin- s hell solution for m = 1 (see IKhourv fc Weltmanl 
(|2004f )) is characterized by the field being stuck at the 
minimum, (f> c , inside the body. We therefore look for so- 
lutions where <f>(0) = 4>i ~ (j) c . The linear approximation 
V c s.<t> = rn 2 ((j> — <p c ) is now valid close to r = with the 
solution 



sinh(m c r) 



m c r 



5 = 



(pi - 4> c 



(47) 
(48) 



We assume that this solution is valid all the way to r — 
R c . For this to be true, the linear term in the Taylor 
expansion of V e s.,p must dominate over the higher order 
terms inside the body. Since (f> is increasing in < r < R c 
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the largest value of 
to the condition 



,\ occurs at r = R c , and leads 
« iZ i 01 ( 49 ) 



Outside the body the linear approximation V e g * 



bb) is valid with the solution 



(50) 



where we have assumed m^Rc < 1 as would be the case 
in most interesting cases. Matching the two solutions at 
r = R c gives us 



.4 



1 



tanh(m c i? c 



TOv.iL 



4> c cosh(m c -R c ) 
The condition Eq. becomes 
4>b - 4>c ( tanh(m e iT > c ) 



(51) 
(52) 

(53) 



For i > 1 we have tanh(x) ~ 1, and by using </>£, <f> c 
we find that the condition above is satisfied for all 



m c R c > — 

<Pc 

The far-away field can now be written 

, , Qcs Mi _ m 



where 



QoS 
£th 



4ttM p1 r 

= 3/3, < p b Mpieth 
4>b - 0c 



6/3,, 



(54) 

(55) 

(56) 
(57) 



and $ c is the gravitational potential for the body. The 
thin-shell factor e t h is on the same form as found in 
iKhourv fc Weltmanl (|2004l ), but it does not have the ge- 
ometrical interpretation as an explicit thin-shell. It is 
however this factor which determines the suppression of 
the fifth-force. 

Comparing the effective coupling Eq. (for 'large' 
bodies: e t h "C 1) with the corresponding expression 
Eq. (|45f (for 'small' bodies: e t h 3> 1) we see that 



{cfi 



= 3e th < 1 



(58) 



This shows that the chameleon force, relative to gravity, 
between two bodies is suppressed as long as one (or both) 
of the bodies have a thin-shell, and demonstrates that the 
chameleon mechanism is present in this model. 

4.1.2. Lunar Laser Ranging 

We will restrict our attention to tests of the equiva- 
lenc e principle using Lunar Laser Ranging (LLR), see 
e.g. [Will (1993). LLR measures the free-fall acceleration 
of the moon and the earth relative to the sun. The ac- 
celeration induced by a fifth force with the field profile 



4>(r) and effective coupling Q oS is amh = |<9effV0|/M p i. 
In most interesting cases, m^ 1 > lAu, the chameleon is 
a free field in the solar-system. This leads to the follow- 
ing constraint (see lKhoury fc Weltmanl (|2004l ) for a more 
detailed derivation) 



2|On 



2Q c Q ff |Q- ff -Q® | < 10 



-13 



(59) 



where Q%, Q^Jg and Qf s is the effective coupling for 
the sun, moon and earth respectively. The background 
density in the solar-system is pb ~ 10 -24 g/cm 3 corre- 
sponding to the average matter (dark and cold) density 
in our galaxy. 

The resulting bounds for m = 2 and m = 3 are shown 
in Fig. ©. 

The strongest constraints on chameleon models (with 
natural para meters) typica l ly com es from the Eot-Wash 
experiment (jKapner et al.l 120071 ). As we can see in 
Fig. ([5]), the interesting part of the parameter space is 
n <C 1. In this regime the test- masses used in the Eot- 
Wash experiment will typically not have thin-shells and 
the experiment cannot, by design, detect the chameleon. 
Nevertheless, the LLR constraints are good enough to 
constraint this part of the parameter space where inter- 
esting cosmological signatures take place. 

4.2. Cosmological constraints 

Due to the conformal coupling, Eq. ([2]), of ^ to matter, 
a constant mass scale mo in the Jordan-frame is related 
to a 0-dependent mass scale m(</>) in Einstein- frame by 
m((f>) = moe^^' . A variation in <fi leads to a variation in 
the various masses 



A/3(0) 



(60) 



BBN constrains the variation in m(</>) from the time of 
nucleosynthesis until today to be less than around 10%. 

If the field has settled at the minimum before BBN, the 
resulting bound turns into /3(</> to day) ^ 0.1 since /3(<j>) is 
an increasing function of time. When the field is not at 



BBN 



< 



the minimum at BBN we must also require 
0.1. 

For our models this constraints only the parameters 
where the coupling satisfies Q^o — lA&oda -^pil ^ 1> m 
which both the background and the matter perturbations 
are completely similar to ACDM. The bounds for Model 
B Eq. jn]| are shown in Fig. (gj). 

Another important restriction on chameleon theories 
comes out from considering the isotropy of the Cosmo- 
logica l Microwave Background (CMB) (Hinshaw et al. 
I2003f l. A difference in the value of cf> today and the value 
it had during the epoch of recombination would mean 
that the electron mass at that epoch differed from its 
present value by -^p^ ~ A/3(</>). Such a change in m e 
would, in turn, alter the redshift at which recombination 



occurred, z roc : 



A.: 



-rc 



A/3(0) 



(61) 



WMAP bounds z roc to be within 5% (at 2a, 23% at 4cr) 
of the value that has been calculated using the present 
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Fig. 4. — Local gravity constraints for the Model B Eq. 11141 
when n = 6 (top) and n = 10 (bottom). The horizontal line 
shows Af = Moe = 10 _3 eV. The vertical line shows A = 1. 
The dashed line shows Q^o = \fl><f>t a a ay l-^pl = 1- The labels 
'Irvine' and 'Eot-Wash' refers to test of the Newtonian gravita- 
tional law, 'Casimir' refers to Casimir experiments, 'PPN' refers to 
test of post-Newtonian gravity in the solar-system and 'BBN' to 
constraints from Big Bang Nucleosynthesis. 



day value of m e (|Nagata et al.ll20 04t). Denoting 0o, <?W 
and </>bbn with the field value today, at recombination 
and BBN respectively. Then /3(0o) > /3(<?W) > /3(0bbn) 
and the CMB bound is weaker then the bound coming 
from BBN. 

5. PERTURBATIONS 

In this section we will study the growth of perturba- 
tions in chameleon models. We start by consider the 
general scalar-tensor model given by the action Eq. ([1]) 
with a universal matter coupling-function P{<f>) and po- 
tential V{4>). In deriving the perturbation we will work 
in units of M p \ = -y==? = 1- For simplicity, we will 
consider the Jordan-frame matter-density satisfying 

pm + (3H - 0,4 </>) p m = (62) 

since this choice will simplify the field equation. In terms 
of the Einstein-frame density p m F this choice corresponds 
to p m = e^^p m F . This is just a matter of convenience 
since ss 1 in the late universe whenever the the- 

ory satisfies the BBN bounds. With this choice the field 
equation reads 

4> + 3H<j> + V,+ +13,4, p m = (63) 

The most general metric in a FLRW space time with 
scalar perturbations is given by 

ds 2 = -(1 + 2a)dt 2 - 2aB. i dtdx i (64) 
+a 2 ((1 + 2t/0<% + 2j Aj ) dx l dx J 

where the covariant derivative is given in terms of the 
three-space metric which in the case of a flat background 
reduces to Sij. We decompose the field cj> into the back- 
ground and perturbations parts: </)(x, t) = <fi(t)+5(j)(x., t). 
The EM-tensors of non-relativistic matter can be decom- 
posed as 

^0 = -Pm(l + 5 m ), if = -p m V,i (65) 

where v is the peculiar velocity of non-relativistic matter 
and 5 m is the matter-density perturbations defined by 

in Pm _ 5p 

o m = v — m the co-moving gauge (6o) 

Pm Pm Pm 

The equation determining the evolution of the per- 
turbations follows from the Einstein-equ ations. In 
the gauge- ready formulation (jHwand 1 199 ID , the scalar 
perturbations equations are 
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3H8(j) + (V44 - -^)^4> + &i44 PmScj) + 2aV,4, + (3,4 (2ap m + 5p m ) - <j>{a - 3Ha + k) 



Sp m + 3H5p m - p rn Ik - 3Ha + ~^y\ - /?,</. {p m S<( + Sp m 4>) - (3,^ p 

A 1 / 

k + 2Hn + 3aH + —a - - I Sp m - 4a</> 2 + 4(f>5(j) - 2V^§ 



A 1 / ... 

A 3 

K+ -2%- -(p m V + 4>64>) 

a 2 

v - a + (3,4 (4>v - 54>) 
X + H x — a — tp 



(67) 
(68) 

(69) 
(70) 

(71) 

(72) 
(73) 



with 

X = a(B + ay) (74) 

K = 3(-ijj + Ha)- (75) 

and A being the co-moving covariant three-space Lapla- 
cian. In the list of equations above Eq. (|6"T|) is the scalar 
field equation of motion, Eq. (|68[) the continuity equa- 
tion, Eq. (|69f the Raychauhuri equation, Eq. (|70|) the 
ADM energy constraint, Eq. (f7"Tjl the momentum con- 
servation constraint and Eq. (f73|) the ADM propagation 



equation. In these equations we have not yet fixed the 
gauge-degrees of freedom. The choice of a gauge will 
simplify the system and we will work in the so-called co- 
moving gauge (v = 0). This gauge leaves no residual 
gauge freedom and we can solve the system for the two 
variables {5<p, 8 m ) directly. 

From Eq. ([71]) we have a = -(3,4 ( 50. So lving Eq. (j6"9"|) 
for k and inserting this into Eq. (|67l69p we find, af- 
ter transforming to Fourier space, the following equations 



S m +2HS n 



2 PmUn 



- (3,4 [6H 2 + 6H- — + 2<j> 2 } 



([3,44 [2H<p - Vegrf] + /3, 4,4,4, <f> 2 ^j - f},4> S4> - 5<t> (5(3,4 H + 2<p + 2/3,44 ^) 



= 



(3H + 2/3,, 



+ P,4> p m S m - #m + ( m| + - 2 P'4 v eff,4 + 2(3 



(76) 
(77) 



where k is a co-moving wavenumber. 

When the field is slow rolling along the minimum we 
can neglect all terms proportional to and the oscil- 
lating term V e g t< p. The perturbations in <f> will evolve 
more slowly than the perturbations in 5 m for scales deep 
inside the Hubble radius, thus, the term p m /3,4 5 m and 

(m? + %?)84> will dominate over the 5(f> time derivatives 
in Eq. (|77p . Using these approximations, we can simplify 
Eq. (J7BJ| and Eq. (J77J) to 



2HS m = 
-3/3,4,a 



3n tt'2 
2 i L m L1 



1 



(78) 



1+- 



Note that the perturbations in 

n 2 



satisfies 

T2 



< 5m (in 

M p i = 1 units) as long as m'^ » H 2 and /3,^ < 0(1) 
justifying dropping the (^-derivatives in Eq. ((771) . 



Restoring M pl 2 = 87rG and defining Q$ — \(3,^M p i\ 
we can write this first equation on the same form as in 
ACDM 



2H6 m = AirGeffPmSn 



where 



a 



elf 



G 1 



2QI 



(79) 



(80) 



1 + 



The quantity G e ff is seen to encode the modification of 
gravity due to the chameleon in the weak-field regime. 

The chameleon will also exhibit an oscillating term, 
but this term is time-decreasing and hence negligible for 
small redshifts. In some /(i?)-models however, this oscil- 
lating term can grow to infinity because the mass of the 
scalaron is not bounded above. The divergence of this 
mass can be removed by adding a UV-term as shown in 
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iThongkool eTaLl (f2009h . 

In studying perturbations, it is convenient to introduce 
the growth-factor / = rf J°^^ ) ) . In ACDM / 1 at high 

redshifts and / — > 1 in an Einstein-de Sitter universe. It 
is important to find a characteristics in the perturbations 
that can discriminate between different DE models and 
the ACDM. Writing the growth factor as 



(81) 



can be a parametrization that is useful for this purpose, 
see e.g. IGannouii et all (|2009allbh . In ACDM we have to 
a good accuracy 7 « 0.55 for low redshifts. Of course in 
some models 7 will vary to much for it to be considered 
a constant, and there can also be a scale dependence, so 
we should write 7 = 7(2:, fc). 

We will be most inter ested in scales k relev ant to the 
galaxy power spectrum (|Tsuiikawa et al1l2009f) 



O.OlhMpc- 1 <k< 0.2hMpc 



(82) 



where h = 0.72 ± 0.08 corresponds to the uncertainty 
in the Hubble factor today. These scales are also in the 
linear regime of perturbations. 

It is also convenient to introduce the length scale of 
the perturbations A p = and the length scale of the 
chameleon A^ = . In Eq. ([80]) we have two asymptotic 
regimes: 



Gcff 



G 

G(l 



A p A,; 
Ap < A,; 



(83) 



In the GR regime, X p 3> A^, the perturbations show no 
deviation from ACDM. 

In the scalar regime, X p <C A^, however, the matter- 
perturbations will feel a stronger gravitational constant 
than in GR. The coupling is in general a dynam- 
ical quantity, which will increase with time when the 
chameleon follows the minimum. Thus, when we reach a 
time where Q$ > 1 the perturbations will start to grow 
with increasing amplitude and will quickly enter the non- 
linear regime . 

See iBaldl (l2010h: ILi et all (I2010alfbh: ILi fe B arrow 
(|2010MLi fe Zhaol(|2009f k lZhao et al.l (|2011l )~ TB"aldi et al. 
(|2010f ) for a numerical N-bo dy simulation a nalysi s 
with this t y pe of mo dels, and | Nimes fc Motal IpOOl : 
iMota et al.l (|2008alTbT ); iBrax et al.l (|2010f ) for a study of 
the spherical collapse in cosmological models with a time 
dependent coupling between dark energy, dark matter 
and other matter fields. 

5.1. The critical length scale A^ 

In order to study the perturbations more closely, we 
look at the value of the critical length scale today A^.o 
for our models. 

The critical length scale for Model A Eq. (|13jl satisfies 



A 



0,0 



10 



(" + 2) 
2 <™ + l)pC 



(84) 



where Q^q = |/3,^ toda M p i\. In order for this length scale 
to affect the matter perturbations: A^o = O(lMpc) and 
Qt/,0 ~ 1, we need to impose n < 0.5. See Fig. [S] for a 
plot of the growth factor j(z — 0). The plot shows the 
two regimes: 



• (i): Phase space where 7 < 0.50 for all relevant 
scales. The perturbations are in the scalar regime. 

• (ii): Phase space where 7 sa 0.55 for all relevant 
scales. This is the GR regime. 

For our Model B Eq. ([T4"|) . the critical length scale is 
given by 



(n- 2) 

M p \ \ 2 <™- 1 ) 2 ( "~ 2 i 27T 
n(n + m)\^J~) Q ^"~'m (85) 



which gives 



A 



4>,o 



10" 



(86) 



In order to have A^o = O(lMpc) together with a cou- 
pling Q^o of the order of unity we need n < 0.5. However, 
for n < 2 the model is no longer a chameleon according 
to our definition in section 1: the range of the field is 
shorter in the low density cosmological background than 
in a high density environment. This also means that local 
gravity bounds will most certainly be violated. 

The only way to increase A^o up to a mega-parsec 
value is by decreasing the coupling strength. For n > 2 
we need 



Qm < 10 



-10 



(87) 

which is to small to significantly affect the growth of the 
perturbations. Thus the perturbations in Model B are 
always in the GR regime. 

It is only in Model A Eq. (fT3|) that we can have in- 
teresting signatures on the matter perturbations. But 
after imposing local gravity constraints we find that the 
perturbations are confined to be in the GR regime with 
no signature on the matter perturbations or on the back- 
ground evoluti on relative to ACDM. T his agrees with the 
result found in IGannouii et all (|2010f ) (for m = 1). The 
only way to have observable signatures in these models is 
to restrict the coupling to dark matter only, and thereby 
avoiding the local constraints. 

The different regimes shown in Fig. (JSJ) have been de- 
rived by considering an universal coupling. Since dark 
matter is dominating over baryonic matter at large scales 
the regimes in this figure is expected to be similar if we 
restr ict the coupling to dar k matter only. 

In IGannouii et all (|2010T ) . the chameleon model 



HW M p i 

y(0) = M 4 (l-/i(l-e _1 tr) n ) 



(89) 



where < n < 1 and < n < 1, was found to 
have observable signatures on the growth of the mat- 
ter perturbations even when local gravity constraints 
was taken into account. We note that this potential 
do not directly generalize to a more general coupling- 
function f3(4>) = (107) ■ This is because the require- 
ment < n < 1, which is required in order to have a 
positive definite mass of the field, leads to violation of 
local gravity constraints for parameters where signatures 
are present. 
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Fi g. 5 . — The two regimes for the growth factor 70 in Model A 
Eq. (1 1 3 1) : (i) 7 < 0.50, and (ii) 7 0.55 for all relevant scales. In re- 
gion (i) we have a significant dispersion between scales. Above, the 
quadratic coupling m = 2, and below the cubic coupling m = 3. 
The label 'Local constraints' shows the allowed region from ex- 
perimental tests of the equivalence principle (Lunar Laser Rang- 
ing). Due to numerical difficulties the pe rturbations have been 
integrated using th e approx imation Eq. 1781 1 rather than using the 
full equations Eq. j76t-fT7l 

6. SUMMARY AND CONCLUSIONS 

We have discussed the cosmological evolution of 
chameleon models with power-law couplings and power- 
law potentials. The chameleon follows the attractor so- 
lution <p = <^min as long as ra^ 3> H, and the attractor 
can be reached for a large span of initial conditions. In 
fact, non-constant couplings can allow for a larger off- 
set from the attractor in the early universe and still be 
in agreement with BBN bounds. Along the attractor 
the chameleon is slow rolling and can account for the 
late time acceleration of the universe. The background 



evolution is however found to be very close to that of 
ACDM, as found in many other similar m odels (see e.g. 
iFaulkner et~ai1 (|2007t ); IBrax et all (|2008l )V The mass- 
scale M in the potential is fine-tuned in the same manner 
as a cosmological constant, thus not providing a solution 
to the fine-tuning problem. This is however motivated 
and required by local gravity experiments. 

Even though the background expansion is very close 
ACDM, the growth of the linear perturbations can be 
quite different. The reason is the fifth-force acting on 
both dark matter and baryons, which leads to a different 
growth rate of matter perturbations on cosmic scales. 

For our Model B Eq. (fTl| . the linear perturbations 
are not affected by the chameleon since the field-range 
is in general too small compared to cosmic scales, or, if 
the field-range is large enough then the coupling is too 
small to produce observable effects. Otherwise the model 
would be ruled out by local gravity constraints. 

In Model A Eq. (TT5|) , the range of the chameleon can be 
large enough as to affect the matter perturbations, which 
leads to a growth-rate different from ACDM. Since the 
coupling, in general, varies with time, we will also have a 
dispersion for scales within the linear regime. However, 
for this to be the case, it must be emphasized that lo- 
cal gravity constraints force us to have a gravitational 
coupling of the chameleon field to dark matter only. By 
neglecting the coupling to baryonic matter, the purpose 
of the chameleon mechanism are lost. But with this con- 
sideration, the growth of matter perturbations can in 
principle allow us to discriminate between Model A and 
ACDM. 

Naively, one would expect that introducing a more gen- 
eral coupling would lead to a more richer phenomenology. 
The reason why this is not the case in the models consid- 
ered here is because local constraints are so stringent it 
forces the chameleon to stay at the local minimum almost 
everywhere in space. This again means that the kinetic 
term of the chameleon is suppressed and the dynamics of 
the theory is determined solely by the interplay between 
the potential and the coupling. In the power-law mod- 
els considered here this in terplay is very simi lar to the 
standard chameleon model IBrax et all (|2010b[ ). 

If future galaxy surveys manages to pin down the mat- 
ter power spectrum to a greater accuracy, and detect a 
clear deviation from ACDM (especially if a scale depen- 
dence in the growth index 7 is detected), then it would 
be interesting to see how good these chameleon models 
can fit the data. 

Even though these models, for most parameters 
choices, do not leave an imprint on the linear matter 
perturbations, it may have an impact on the small scale 
structure formation. It would be interesting to investi- 
gate the effect of these models on the non-linear regime of 
structure formation by means of high resolution N-body 
simulations. We leave this for future work. 
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